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Unsteady Duct flows 


Unsteady plane shear flows are experienced 
by fluids where the driving forces are 
suddenly or continuously changed. A fluid 
initially at rest will for example respond to a 
sudden increase in the pressure gradient and 
the volume flow through the duct will then 
increase monotonically until steady state is 
reached. 
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Starting flow in circular pipes 
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Matlab Code 


plot3(r(1:m),t(1:m),u, linewidth’, 2); 
grid on 


> Clear; 

» clc; 

» r=1:0.1:100; 

> t=0:0.1:100; 

> Um=1: 

> JO=1; 

> J1=1; 

ኑ k=40; 

» lamda=2; 

> lamda=lamda*ones(1 ,k); 

» m=min(length(r),length(t)); 
> for ii=1:m 

> u(i)=compute u [ t(r(ii),t(ii),Um JO J1,lamda); 
> end 

» figure; 

2 

2 


» function u=compute u r t(r,t,Um,JO,J1,lamda) 
» u=Um*(1-r.A2); 

» tmp-28*Um*JO*lamda'r; 

» tmp=tmp./(lamda.A3); 

» tmp=tmp/]1; 

» tmp=tmp./lamda; 

» tmp-tmp.*exp((-lamda.^2).*t); 

< sumTmp=sum(tmp); 

» u=u-sumImp; 

» end 


» function u=compute u r t(r,t,Um,JO,J1,lamda) 
» u=Um*(1-r.A2); 

» tmp-28*Um*JO*lamda'r; 

» tmp=tmp./(lamda.A3); 

» tmp=tmp/]1; 

» tmp=tmp./lamda; 

» tmp-tmp.*exp((-lamda.^2).*t); 

< sumTmp=sum(tmp); 

» u=u-sumImp; 

» end 
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WE NEED TO KNOW HOW WE CAN OBTAIN 
THIS SOLUTION IN DETAIL 
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Figure 6.31. Velocity contours for steady unidirectional flow in tubes of rectangular 
cross sechion with width-to-height ratio equal to 1, 2 and 4. 


